Electrostatic screening in fullerene molecules. 
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The screening properties of fullerene molecules are described by means of a con- 
tinuum model which uses the electronic wavefunctions of planar graphite as a starting 
point. The long distance behavior of the system gives rise to a renormalizable the- 
ory, which flows towards a non trivial fixed point. Its existence implies an anomalous 
dielectric constant. The screening properties are neither metallic nor insulating. Al- 
ternatively, the intramolecular screening is obtained from a simple approximation to 
the electronic wavefunctions. Intermolecular effects are also calculated, As a consis- 
tency check, it is shown that the observed polarizability of Cgo is well eproduced. 
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The localization of the electrons within the Cgo spheres implies that electrostatic inter- 
actions cannot be neglected in fullerene crystals. The simplest estimate for the charging 
energy of a C 60 molecule, e 2 /R, where R is the radius, gives a rather large value, ~ 4eV. 
The intermediate size of the molecule puts it in a regime different from an isolated atom 
or a bulk system. Different approaches have been used in the study of electrostatic effects 
0Jiy3,@]. It has also been proposed that they can be the origin of superconductivity in doped 
systems 

In the present work, we use a simple model for the long wavelength electronic properties 
of fullerene molecules [7|||, to analyze the screening properties. It has been shown that it 
gives a reasonable approximation to the electronic levels of Cgo- Furthermore, the scheme 
is sufficiently general to describe other systems with similar topology, like nanotubes and 
fullerene 'onions'. As shown in a previous work ||, the model can be used to study other 
features, like the electron-phonon interaction' with reasonable accuracy. Calculations based 
on this model || for the electrostatic interactions between the highest occupied orbitals (in 
doped systems) are in good agreement with more detailed calculations ||. 

The scheme is based on the fact that the low lying electronic states of an isolated graphitic 
sheet are well approximated by an effective Dirac equation in (2+1) dimensions, using 
Hiickel's theory of conjugated carbon compounds. This description plays the same role 
as the effective mass theories used in the long wavelength regime of other materials. The 
existence of rings with an odd number of atoms in Cgo leads to additional complications, 
which require the introduction of a fictitious gauge field, to account for the frustration that 
these rings induce in the electronic states. 

We use the preceding scheme to analyze first the screening properties of an isolated 
graphite plane. The fact that the dielectric constant is anomalous can be inferred from the 
semimetallic nature of the system. The density of states at the Fermi level is zero, and there 
is no Thomas Fermi screening at long distances. On the other hand, there is no energy gap, 
so that the system does not behave like an insulator. We can write an effective action for 
the system, coupling the Dirac electrons to the electromagnetic field in the standard way: 

S = J d 2 rdt \P"(7o<9 + v F j • a )* ~ i ~ c J d 2rdt *(7oA + • A )* (1) 

The action is not Lorentz invariant, as the velocity of propagation of the fermions, vp, is 
different (and much smaller) from that of the photons, c. Moreover, note that the Fermi 
fields are defined in a two dimensional plane, while the vector potential exists throughout the 
entire three dimensional space. In this way we ensure the correct interaction at long distances 
between charges in the plane, e 2 /|r — r'| . A consistent coupling of the electromagnetic field 
is attained, in the quantum theory, by performing the gauge fixing so that the spatial 
component transverse to the plane decouples from the electronic system. This is achieved 
precisely by the Feynman gauge, in which has the same direction as the current and each 
component does not mix with the rest 

r d 3 k du P *k-(r-r')-Mi-0 
{ T M t. r)M<, r')> - ~i V / ^ P) 

From dimensional counting, it is easy to conclude that the coupling constant e is di- 
mensionless (in units in which h — c — 1), so that perturbation theory should give rise 
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to logarithmic corrections, which can be included in renormalized parameters in the stan- 
dard way. We are interested in the low energy properties, which are given by the infrared 
fixed point of the theory. Our initial conditions are fixed by the known bare parameters, 
vp and e at the upper cutoff where our continuum model is valid. This cutoff, in energies, 
is determined by the bandwidth of the tt electron band that we are considering, which is 
proportional, within Hiickel theory, to the interatomic hopping matrix element of the hamil- 
tonian, t ~ 2.2eV [PH|]. In space, the relevant cutoff is the interatomic distance, a = 1.4A. 
In terms of these parameters, Vp = 3ta/(2h). 

Regularizing the theory by analytic continuation to dimension d = 3—e |TT|, the divergent 



contributions to the quantum effective action are given, to lowest order in e 2 /hc, by the 
diagrams in Figure 1. It can be shown that the two dimensional matter fields do not 
renormalize the scale of the electromagnetic field. The self-energy contribution from Figure 
la gives rise to the renormalization of Vp, Vp = Z Vp (vp)n, as well as to the wavefunction 
renormalization, \& = r. In terms of hypergeometic functions, we obtain: 

Z, y I — — { n—F ! -. -: ~ + i - 2tt^ [ 1 - 
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The correction to the vertex from Figure lb gives the renormalization of the electric charge, 
e = Z e eR, and provides at the same time an alternative way of computing the renormalization 
of Vp. The value of Z VF from the vertex correction coincides precisely with the above 
expression from the self-energy diagram, ensuring the renormalizability of the model. The 
second remarkable point is that, to this perturbative order, Z e turns out to be equal to 1. 
This means that both pieces of the action (1) are affected by the same global renormalizatiton 
and that, therefore, the gauge invariance of the theory is maintained at the quantum level. 

In the relevant physical limit, vp/c — > 0, the structure of (3) is greatly simplified, and 
we obtain: 

Z Vp = 




Z* = l + — (4) 



Note that the expansion, which originally was on the fine structure constant, e 2 /(hc) 
is changed into an expansion on e 2 /{Tivp) cx e 2 /{ta). This was to be expected on physical 
grounds, as in the limit when the propagation of the electromagnetic field is nearly instan- 
taneous, this is the only measure of the relative strength of the electrostatic and kinetic 
energies. This parameter is greater than unity, although the convergence of the series is 



3 



improved by the presence of dimensionless constants in the denominators of all terms. The 
fact that the electric charge is not renormalized is true to all orders, in the limit vpjc — > 0. 
The only coupling which remains is longitudinal, and the matrix elements for the relevant 
transitions tend to zero for small momentum transfers. 



Applying standard renormalization group methods [11|, we obtain the variation of the 



renormalized coupling ta/e as a function of the energy scale \x 




The computation of higher perturbative orders simply amounts to adding higher powers 
of e 2 / (ta) to the right hand side of the equation, and does not alter the qualitative behaviour 
of the renormalized coupling. From (5) it is clear, in fact, that perturbation theory becomes 
more accurate as the spatial scale increases in the infrared regime. 

The scaling of vf changes the density of states al low energies, which behaves as: 

N(e) ~ | e |i+e 2 /(4M+o(e 4 /(M 2 ) ( 6 ) 

The scaling outlined above is limited, in a fullerene sphere, by its radius, which acts 
as an infrared cutoff. In a doped compound, another low energy scale arises: the Fermi 
energy associated with the extra electrons, or a length oc kp l . In the cases of interest, the 
filling factor is small, 3/60, so that the relevant cutoff remains the radius of the molecule 
(hp 1 ~ 20a). Note also that, in doped systems, the existence of a finite density of states at 
the Fermi level allows us to define a Fermi-Thomas screening length ( k FT ~ e 2 kp/{taR) ). 
Again, for the relevant values of hp, this length is greater than the radius. Hence, metallic 
screening plays no role on the intramolecular screening in doped Ceo- 

The absence of renormalization of e 2 , mentioned earlier, implies that the charging energies 
scale as e 2 /R, although with a proportionality term reduced by the finite corrections. 

We can also use our continuum aprroximation to the electronic structure of the fullerenes 
to study their electrostatic properties in a different way to the one discussed before. In prin- 
ciple, the knowledge of all the (discrete) electronic levels, allows us to use standard second 
order perturbation theory, and calculate the polarization of the inner levels as fucntion of 
the extra electrons added to the molecule. In this way we can estimate the attractive in- 
teraction between two additional electrons, induced by the polarization of the other levels. 
The level spacings scale like ta/R, so that, on dimensional grounds, this interaction behaves 
as e A /(taR), while the direct Coulomb repulsion is oc e 2 /R. Thus, the balance between the 
two depends crucially on the dimensionless prefactor. For the wavefunctions of the lowest 
unoccupied triplet, this value is given by a convergent sum, which involves transitions from 
multiplets which are higher and higher in energy. The convergence of this sum has the same 
origin as the finiteness of the diagrammatic corrections to e 2 in the previous analysis. From 
the computation of this series, we find a value of 0.14. Hence, the polarization of the inner 
levels in C 6 o reduces by a significant factor ~ 50 % the bare Coulomb interactions. We can 
also compute the effective interaction between two additional electrons in the lowest triplet. 
This interaction is attractive, and reduces strongly, but does not cancel, the bare repulsion. 

Finally, we can compute the interaction between electrons located in different molecules. 
For two large molecules at close distance, the most convinient starting point is tha analysis 
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of two planes of graphite. The screening of the electrostatic interactions is described by 
the scaling given in eq. (5), until a scale proportional to the interplane distance is reached. 
Then, additional screening is induced by the mutual polarization effects of the two planes. 
This effects leads to an interaction, which, to lowest order, goes like e 4 x A/(taD 3 ), where 
A is the area of the planes, and D the interplane distance. 

At sufficiently large distances, each molecule is described by its polarizabilty in a standard 
way. It can be computed from the knowledge of the electronic transition energies, and the 
corresponding matrix elements. For a neutral, spherical fullerene, we find: 



e 2 R 3 
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(7) 



Inserting the appropiate parameters for Cgo, we obtain P = 90A, which is in very good 
agreement with the experimental value, 80 A [|lt|l2|, as well as with more sophisticated 
theoretical calculations P, |T3| , |l4|| . 

From the scaling of the density of states (6), we can also infer the role of intermolecular 
tunneling. The anomalous terms in the exponent reduce the effective tunneling at low 
energies, in a similar way to the case of a ID Luttinger liquid [15-17]. If the correction 



to the non interacting exponent is greater than one, it can be shown that intermolecular 
tunneling becomes an irrelevant perturbation at low energies. In the present situation, that 
is not tha case. However, the surviving effective tunneling should be greatly reduced with 



respect to its bare value fl7| : 
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(8) 



In summary, we have studied the screening properties of fullerene molecules. The perfect 
semimetallic nature of isolated graphite planes leads to unusual features, like the absence 
of metallic screening, although the system does not behave like an insulator. There is a 
significant reduction of the intramolecular charging energies, and an associated reduction in 
the effective intermolecular hopping. We have also shown how to compute the most relevant 
intermolecular electrostatic effects. 

From a fundamental point of view, we have found that the low energy properties of the 
system give rise to a non trivial, renormalizable, field theory, which looks like a simplified 
version of QED. The anomalous dimensions acquired by the propagators resemble, in many 
ways, the one dimensional Luttinger liquids. 
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FIGURES 

FIG. 1. Low order diagrams computed in the text, a) Self energy correction, b) Vertex 
correction. The convention used for the momenta and energies is that of the limit vpjc — > 0. 
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